Universal critical behavior of aperiodic ferromagnetic models 
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We investigate the effects of geometric fluctuations, associated with aperiodic exchange inter- 
actions, on the critical behavior of q-state ferromagnetic Potts models on generalized diamond 
hierarchical lattices. For layered exchange interactions according to some two-letter substitutional 
sequences, and irrelevant geometric fluctuations, the exact recursion relations in parameter space 
display a non-trivial diagonal fixed point that governs the universal critical behavior. For relevant 
fluctuations, this fixed point becomes fully unstable, and we show the apperance of a two-cycle which 
is associated with a novel critical behavior. We use scaling arguments to calculate the critical expo- 
nent a of the specific heat, which turns out to be different from the value for the uniform case. We 
check the scaling predictions by a direct numerical analysis of the singularity of the thermodynamic 
free-energy. The agreement between scaling and direct calculations is excellent for stronger singu- 
larities (large values of q). The critical exponents do not depend on the strengths of the exchange 
interactions. 



Quenched disorder may change the critical behavior of 
ferromagnetic spin systems |l||. Although disorder pro- 
vides the more usual examples, there are alternative ways 
to break translational invariance. For instance, the inter- 
est in the study of quasicrystals inspired a number of 
proposals of deterministic schemes to build spatially ape- 
riodic structures. In analogy with the Harris criterion to 
account for the effects of quenched disorder. Luck Q de- 
veloped a heuristic reasoning to gauge the relevance of 
geometric fluctuations (associated with aperiodic inter- 
actions) on the critical behavior of ferromagnetic models 
on Bravais lattices. In a recent publication Q, Luck's 
criterion has been exactly derived for a g-state ferromag- 
netic Potts model on a diamond-type hierarchical lattice. 
Now we revisit this problem to show the existence of an 
attractor that gives rise to a novel class of critical behav- 
ior. 

Many investigations of aperiodic classical pi and quan- 
tum H Ising chains, and two-dimensional Ising models 
with aperiodic layered interactions H , use the formalism 
of substitution rules on alphabets for the construction of 
infinite aperiodic sequences of elements, which arc then 
associated with coupling constants along a direction of 
the lattice. We have taken advantage of the structure of 
hierarchical lattices |q| to build layered aperiodic mod- 
els pi which are then amenable to simple (and exact) 
renormalization-group calculations. In parameter space, 
there is always a "diagonal" , non-trivial, fixed point of 
the renormalization-group recursion relations. For rele- 
vant geometric fluctuations, this diagonal fixed point, as- 
sociated with the critical behavior of the uniform (equal 
couplings) model, becomes fully unstable, and therefore 
cannot be reached from any nonuniform initial conditions 
in parameter space [Q . The critical behavior of the ape- 
riodically perturbed system, in case of a phase transition, 
should then be governed by another attractor in param- 
eter space. 

Relevant aperiodicity is believed to drastically weaken. 



or eventually suppress, the critical singularities p]. Re- 
cently, convincing Monte Carlo evidence |10| has been 
presented to indicate that relevant layered aperiodicity in 
the 8-state Potts model on the square lattice drives the 
phase transition from first to second-order (with critical 
exponents independent of the strength of the couplings). 
In the present work, we also detect a weakening of the 
critical singularities for relevant geometric fluctuations. 
Besides being numerically exact, our results can be re- 
garded as equivalent to a Migdal-Kadanoff approxima- 
tion for the analogous problems on Bravais lattices, and 
may thus be put in perspective with the Monte Carlo 
findings |lO|| . 

The g-state Potts ferromagnet is defined by the Hamil- 
tonian 






(1) 



where Ci = 1, 2, . . . , g for all sites of a lattice, Jij > 0, 
and the sum over {i,j) refers to nearest-neighbor pairs 
of sites. We assume that the couplings can take only 
two values, J a and Jg, associated with the sequence 
of letters A and B of an aperiodic substitutional word. 
To generate this sequence, we can use, for example, the 
successive iterations of a period-doubling rule, given by 
{A, B) — > {AB^ AA). In FigjI], we show a simple diamond 
hierarchical lattice with aperiodic interactions according 
to this sequence (the letters, and the corresponding cou- 
pling constants, are chosen to mimic a layered struc- 
ture). In general, we may consider basic "diamonds" 
with m branches and b bonds along each branch, and 
generate hierarchical lattices with intrinsic, or fractal, di- 
mension, D ~ ln(rn6)/ ln(5). In each one of these struc- 
tures, aperiodicity may be implemented by a substitu- 
tion rule of the form {A,B) -> (A"iB''-"i, ^"^B^-"^), 
with Q < ni,n2 < b. These sequences are characterized 
by a substitution matrix with eigenvalues Ai = 6 and 
A2 = ni — 712. The asymptotic form of the fluctuations 



in the number of letters depends on the wandering expo- 
nent, 



log 1^2 I 
logAi 



(2) 




FIG. 1. Three successive generations of the simple diamond 
lattice {b — m = 2). 

Now we decimate the internal degrees of freedom of 
the diamonds to write exact recursion relations for the 
reduced couplings. There is always a non-trivial fixed 
point on the diagonal of parameter space, which is asso- 
ciated with the critical behavior of the uniform model. 
This fixed point becomes fully unstable for 



w > 1 



D 



(3) 



where the critical exponent a„, associated with the spe- 
cific heat of the underlying uniform model, depends on 
q, m and h (note that the transition is always continuous 
for the ferromagnetic Potts model on these hierarchical 
lattices pfl). Therefore, Eq. (0) is an (exact) alter- 
native statement of Luck's original heuristic criterion of 
relevance of the geometric fluctuations. Irrelevancy of 
the fluctuations corresponds to accessibility of the non- 
trivial symmetric flxed-point, whatever the initial condi- 
tions. In the irrelevant cases, the critical behavior of the 
aperiodic systems remains unchanged with respect to the 
corresponding uniform models. 

The problem now consists in the characterization of the 
new critical behavior under relevant geometrical fluctua- 
tions. Andrade ||l2| has recently used a mapping for the 
succession of partition functions of aperiodic Ising models 
{q = 2) on hierarchical lattices to investigate this ques- 
tion. For irrelevant geometric fluctuations, the critical 
behavior seems to remain unchanged with respect to the 
uniform system. Although the results are not conclusive, 
the critical behavior does seem to change for an Ising 
model on a diamond-type lattice with m = h — Z^ and 
layered aperiodic exchange interactions according to the 
rule {A,B) -^ {ABB, AAA). To go beyond this work, 
we analyzed some families of aperiodic g-state Potts fer- 
romagnets on diamond-like lattices to show that they do 
present a novel critical behavior (different from the uni- 
form case, but universal) under relevant geometric fluctu- 
ations. This novel critical behavior can be exactly traced 



to new features in parameter space, exclusively related 
to the relevant geometric fluctuations. 

Consider the Potts model on a lattice with fe = 2 bonds 
per branch, and interactions according to the period- 
doubling rule. From Eq.(y|), for to = 2, fluctuations are 
relevant if q > 4 + 2 a/2 (for ttt, = 3, the threshold value 
of q turns out to be somewhat lower, and so on). The 
recursion relations are given by 



, , xaxb + q-l 

Xjy 



xa-\- xb + q — 2 



and 



^B 



g-1 



2xA 



(4a) 



(4b) 



where xa,b = exp{qPJA,B), with (3 = l/fc^T. Besides 
the uniform (diagonal) fixed-point, we have found a cycle 
of period two in the relevant regime (that is, the two-cycle 
is present when the non-trivial diagonal fixed point is 
fully unstable). Moreover, the two-cycle appears clearly 
as a bifurcation, its points running away from the sym- 
metric fixed point as q increases. In terms of the second 
iterate of the recursion relations (of which each point of 
the two-cycle is a fixed point), it displays a saddle-point 
character, with stable and unstable manifolds. Suppos- 
ing that this two-cycle is associated with a novel critical 
behavior, we use standard scaling arguments to predict 
the critical exponent a of the specific heat. Let us call 
X the single relevant variable describing the parameter 
space in the neighborhood of the two-cycle. In the ther- 
modynamic limit, we write the reduced free energy per 
bond in the scaling form |0] , 



f{x) = g{x) 



1 



/(^"), 



(5) 



where g{x) is a regular function related to the free energy 
of the decimated spins, x" is a second iterate of the recur- 
sion relations, and h is the linear rescaling factor of the 
renormalization-group transformation (which coincides, 
for diamond- like lattices, with the parameter b of the ba- 
sic diamonds [Q). Note the presence of the factor 6^^, 
related to the need of two iterations to go back to the 
neighborhood of the starting point in parameter space. 
Eq.(H) has the asymptotic solution 



/(^) 
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where x* is one of the points of the two-cycle, Kdc is the 
largest eigenvalue of the linearization of the second iter- 
ate of the recursion relations about any one of the points 
of the cycle, P{z) is an arbitrary function of period one, 
and the critical exponent a, associated with the specific 
heat, is given by 



ln&^ 



2-2 



\ii(mb) 



(7) 



In Acic In Acic 

The values of a predicted by Eq. (M) are unequivocally 
different from the values q;„ for the specific heat exponent 
in the uniform cases, as can readily be seen in Table |, 
for the simple diamond lattice, m — b = 2, with layered 
exchange interactions according to the period-doubling 



sequence. As in the Monte Carlo simulations for the 
aperiodic 8-state Potts model on the square lattice pO[ , 
there is a clear weakening of the critical singularities due 
to the geometric fluctuations (a feature also present in 
disordered systems). Similar conclusions can be drawn 
for fo = 2 diamond-type lattices with different values of 
m and q, and the period-doubling substitution. 



TABLE I. Results for the location of the two-cycle, eigenvalues of the second iterate of the recursion relations about the 
points of the two-cycle and the specific-heat critical exponent, a, as predicted by Eq.M, for some values of q, in the case of 
the m — b = 2 diamond with the period-doubling substitution. The value of the exponent in the uniform {J a ~ Jb) case, a„, 
is also shown for comparison. 
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Eigenvalues of 2nd iterate 
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(5.285. . 


. ,7.642. . 


•) 


(6.697. . 


.,4.750...) 


25 


(6.942. . 


. ,234.34. 


••) 


(39.023. 


..,3.831...) 


100 


(181.71. 


..,5.721. 


■•) 


(13.753... 


.,5151.84... 



3.993. . . 


0.985... 


-0.0022. . . 


0.010, 


4.243. . . 


0.343. . . 


0.0817... 


0.404, 


4.975. . . 


0.074. . . 


0.2720. . . 


0.648, 



We performed the same analysis for the Potts model on 
the b ~ 3 generalized diamond lattice, with interactions 



according to the rule {A, B) ^ {ABB, AAA). Now, the 
recursion relations are given by 



J 



x'a 



and 



XAxl + (q- 1)XA + 2(g - 1)XB +q{q-3) + 2 
+ 2xAXB + iq~ '2)xA + 2{q - 2)xb + q{q - 3) + 3 



x\ + 3{q~l)xA + q{q-3) + 2 



3xj^ + 3{q - 2)XA + q{q-3) + 3 



(8a) 



(8b) 



r 



For all values of m, the condition of relevance, given by 
Eq.(H), is satisfied for g > 2 (including the Ising model, 
q = 2). Again, besides the fully unstable fixed point, we 
have detected the presence of a two-cycle. The same sort 



of scaling analysis has been performed. Table pi gives 
some results for b — 3 and m — 2. We see that the 
weakening of the critical singularities is again indicated 
by these data, the same trend being present for m — 3. 





TABLE IL Same as Table |, for the lattice with b 


= 3 and m - 


= 2, 


with the rule 


{A,B) 


^ {ABB, AAA). 






Location of the two-cycle 


Eigenvalues of 2nd iterate 
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{xa,Xb)i (xa,Xb)2 


Ai 




A2 
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Qu 


2 


(6.446. . . ,135.10. . . ) (34.794. . . ,5.224. . . ) 


3.255. . . 




0.311... 




-1.0358. . . 


-0.902. . . 


7 


(11.469. . . ,1649.66. . . ) (126.59. . . ,8.525. . . ) 


4.755. . . 




0.075. . . 




-0.2981... 


-0.185... 


100 


(57.223. . . ,1121275.02. . . ) (3272.70. . . ,34.683. . . ) 


9.808. . . 




0.001... 




0.4304. . . 


0.898. . . 



The presence of these two-cycles seems to be associ- 
ated with the alternance of two energy scales, given by 
J A and Jb . The iteration of the recursion relations leads 
to the alternative dominance of each one of them (as it 
can already be seen in the behavior of the fluctuations 
in the number of letters A and B along a substitution 
sequence). To check this argument, we performed some 
calculations for the ferromagnetic Potts model on a sim- 



ple diamond lattice, m = b = 2, with aperiodic inter- 
actions according to the four-letter Rudin-Shapiro rule, 
{A,B,C,D) -^ {AC,DC,AB,DB), whose fluctuations 
are known to be relevant even in the Ising case ||9|]. In- 
deed, we have found a cycle of period four, along with 
fully unstable two-cycles, in regions of parameter space 
associated with the symmetries of the sequence. 

To test the validity of the scaling arguments, and of 



the role of the two-cycle as the responsible for the new 
critical behavior, we have performed a direct numerical 
analysis of the singularity of the free energy. In real-space 
renormalization-group calculations, it is well known that 
the (reduced) free energy can be expressed as an infinite 
series |15[. For the Potts model on the 6 = 2 diamond- 
type lattice, with the period-doubling rule, it takes the 
form 



f{xA,XB) = 



1 



^^ (2m)" 



iM^f 



„(») 



q-2 



-ln(2a;^ 



(n) 



q-2 



(9) 



where x^g indicate the n-th iterates of the recursion 
relations, Eq.(|4|). If we assume uniform convergence, 
this series can be differentiated term-by-term, and then 
summed in a computer, to obtain the specific heat per 
bond. We also used direct numerical differentiation as a 
control of this assumption of uniform convergence. The 
critical temperature can be determined with high preci- 
sion by making use of the existence of the trivial param- 
agnetic fixed point at T = 0, corresponding to xa.b = oo, 
which causes the apparent divergence of the series (||) if 
summed without the use of any regularization trick. Fix- 
ing the parameters q and tti, and also the strengths of Ja 
and Js, the critical temperature thus calculated in fact 
places xa,b in the attraction basin of the two-cycle. For 
irrelevant aperiodicity, as well as for the uniform model, 
this method yields a critical temperature that locates the 
system on the stable manifold of the uniform fixed point. 
The singularity in the specific heat can be analyzed 
by a (non-linear) fitting of a function to the data over a 
somewhat arbitrary scaling region. For the uniform and 
irrelevant cases, we have always obtained very good fit- 
tings, in excellent agreement with the values of q;„ pre- 
dicted by the usual scaling theory around the uniform 
fixed point. As it should be anticipated, these fitted val- 
ues did not present any detectable sensitivity on the val- 
ues of J A and Jb- The particular problem of the Ising 
model (g = 2) on the simple m = b = 2 diamond lat- 
tice, with exchange interactions according to a period- 
doubling sequence, had already been studied by Andrade 
||l2| , with the same conclusions. In the relevant cases, the 
situation is much subtler, and demands a more refined 
analysis. For large values of q, in which cases the scal- 
ing analysis predicts positive values of a (although, of 
course, smaller than the corresponding values of q;„), the 
fittings presented excellent agreement with the scaling 
predictions. For weaker singularities (mainly a negative), 
the fitted values were always somewhat bigger than the 
scaling predictions, with better agreement for increasing 
values of q. For m = b = 2, and the period-doubling 
sequence, let us give some results of the fittings of the 
specific heat data to a function of the form A -f i? |i|~", 
where t is the reduced temperature, and the parame- 



ters A and B must not be universal. For q = 7, we ob- 
tained a = —0.005(4), which should be compared with 
a = -0.0022 .... For <7 = 25, we have a = 0.08(2), to 
be compared with the scaling value a = 0.0817. . .. For 
q — 100, we have a ~ 0.27(1), to be compared with 
a = 0.2720 . . .. Even in the cases of disagreement with 
the scaling predictions, the fittings indicate no sensitivity 
on the particular values of J a and Jb, and thus charac- 
terize the universal nature of the critical behavior. For 
weak singularities, the discrepancies in the results cannot 
probably be explained in terms of lack of numerical pre- 
cision, although it is in fact difficult to obtain a fully reli- 
able fitting in these cases. This behavior can be probably 
traced to corrections to scaling that we are not consid- 
ering in the direct application of the scaling ideas to the 
two-cycles. It should be pointed out that the free energy 
takes different values in each point of the two-cycle, in 
such a way that it may be ill-defined in terms of just a 
simple scaling field. 

We have performed the same kind of numerical check 
for a lattice with 6 = 3, and with the aperiodic rule 
{A, B) — > [ABB, AAA). Now, aperiodicity is already rel- 
evant for q — 2 (Ising model), and the two-cycle can be 
found for q > 2, and all values of m. The situation turns 
out to be exactly the same as before. The numerical value 
of the critical temperature indeed locates the system on 
the attraction basin of the two-cycle. The agreement be- 
tween the numerical fittings and the scaling predictions 
for a improves as the critical singularity grows stronger, 
which happens with increasing values of q. In the special 
case q = 2, and tti = 6 = 3, we and Andrade p2| have 
calculated similar values, a = —0.90(9), to be compared 
with the scaling prediction, a = —0.9684 . . .. 

From the numerical calculations, we have observed an 
oscillatory behavior in the specific heat as a function of 
temperature above Tc- The period of these oscillations 
is roughly given by Eq.(|6|), with better agreement for in- 
creasing values of q. These oscillations have also been 
found by Andrade [|l2| , and are well known phenomena 
related to hierarchical structures |lq ]. 

In conclusion, we have given a number of examples of 
ferromagnetic Potts models on diamond-type hierarchical 
lattices to show that irrelevant geometric fluctuations do 
not change the (universal) critical behavior with respect 
to the uniform cases. On the other hand, for two-letter 
substitutional sequences, relevant geometric fluctuations 
give rise to a novel universal critical behavior associated 
with a two-cycle in parameter space. 
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